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Abstract 

A number theoretical model of 1 // noise found in phase locked loops is developed. The dynamics 
of phases and frequencies involved in the nonlinear mixing of oscillators and the low-pass filtering 
is formulated thanks to the rules of the hyperbolic geometry of the half plane. A cornerstone of 
the analysis is the Ramanujan sums expansion of arithmetical functions found in prime number 
theory, and their link to Riemann hypothesis. 
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1. INTRODUCTION 



1// noise, already discovered by Nyquist in resistors at the dawn of electronic age, still 
fascinates experimentalists and theoreticians. This is because none physical principle or 
mathematical function allows to predict easily the observed hyperbolic 1// power spectral 
density. Nevertheless the occurrence of 1// fluctuations in areas as diverse as electronics, 
chemistry, biology, cognition or geology claims for an unifying mathematical principle Q|. 

A newly discovered clue for 1// noise lies in the concept of a phase locked loop (or PLL) 
j^. In essence two interacting oscillators, whatever their origin, attempt to cooperate by 
locking their frequency and their phase. They can do it by exchanging continuously tiny 
amounts of energy, so that both the coupling coefficient and the beat frequency should 
fluctuate in time around their average value. Correlations between amplitude and frequency 
noise were observed 0. 

Fortunately one can gain a good level of understanding of phase locking from quartz 
crystal oscillators used in high frequency synthesizers, ultrastable clocks and communication 
engineering (mobile phones or so). The PLL used in a FM radio receiver is a genuine 
generator of 1// noise. Close to phase locking the level of 1// noise scales approximately 
as 5"^, where a = aK/uB is the ratio between the open loop gain K and the beat frequency 
oJb times a constant coefficient a. The relation above is explained from a simple non linear 
model of the PLL known as Adler's equation 

e{t) + KH{P) sin e{t) = oob, (1) 

where at this stage H{P) = 1, ujb = uj{t) —ljq is the angular frequency shift between the two 
quartz oscillators at the input of the non linear device (a Schottky diodes double balanced 
mixer), and 6{t) is the phase shift of the two oscillators versus time t. Solving and 
differentiating one gets the observed noise level a versus the one cr = Sujb/ujb for the open 
loop case. Thus the model doesn't explain the existence of 1// noise but correctly predicts 
its dependance on the physical parameters of the loop ^. 

2. LOW PASS FILTERING AND l/f NOISE 

Besides one can get detailed knowledge of harmonic conversions in the PLL by accounting 
for the transfer function H{P), where P = ^ is the Laplace operator. If H[P) is a low pass 
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filtering function with cut-off frequency /c, the frequency at the output of the mixer + filter 
stage is such that 

^ = fB{t)/fo = qi\i^ - Pi/qi\ < fc/fo, Pi and qi integers. (2) 

The beat frequency fsit) results from the continued fraction expansion of the input fre- 
quency ratio 

= fif)/k = [ao;ai,a2,---ai,a, •••], (3) 

where the brackets mean expansions ao + l/{ai + l/{a2 + l/ .. . + l/{ai + l/{a + .. .)))). The 
truncation at the integer part a = [-j^] defines the edges of the basin; they are located at 
z/i = [ao; ai, 02, ... , a^, a] and z/2 = [oq; ai, 02, ... — 1, 1, a] 4]. The two expansions in z/i 
and z/2, prior to the last filtering partial quotient a, are the two allowed ones for a rational 
number. The convergents Pi/qi at level i are obtained using the matrix products 
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FIG. 1: The intermodulation spectrum at the output of the mixer+filter set-up 



Using dU) one can get the fractions ui and 1^2 as ui = — and z/2 = so that, with 

the relation relating convergents {piqi-i — Pi-iqi) = (— 1)*^^, the width of the basin of index 
i is I z/i — z/2 1 = —7 — -, ^ N — ^— whenever a > 1 (see also Fig. 1). 

In our previous publications we proposed a phenomenological model for 1// noise in the 
:'LL, based on an arithmetical function which is a logarithmic coding for prime numbers 
2|, ^. If one accepts a coupling coefficient evolving discontinuously versus the time n as 
K = KoA{n), with A(n) the Mangoldt function which is ln(p) if n is the power of a prime 
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number p and otherwise, then the average couphng coefficient is Kq and there is an 
arithmetical fluctuation e{t) 



te{t) = - ln(27r) - i ln(l - r^) - ^ 



p p ' 



(5) 



The three terms at the right hand side of te{t) come from the singularities of the Riemann 
zeta function Ci^)-, that are the pole at s = 1, the trivial zeros at s = —2/, / integer, and 
the zeros on the critical line 3fJ(s) = ^Q. Also the power spectral density roughly shows 
a 1// dependance versus the Fourier frequency /. This is an unexpected relation between 
Riemann zeros (the unsolved Riemann hypothesis is that all zeros should lie on the critical 
line) and 1// noise. 

We improved the model by replacing the Mangoldt function by its modified form h{n) = 
A{n)(j){n)/n, with (f){n) the Euler (totient) function ^. This seemingly insignificant change 
was introduced by Hardy in the context of Ramanujan sums for the Goldbach conjecture 
and resurrected recently by Gadiyar and Padma for analyzing the distribution of prime pairs 
0]. Then by defining the error term eB{t) from 



Bit) = J2Kn)=til + eB{t)), 



(6) 



n=l 



its power spectral density Ssif) — jk^ exhibits a slope close to twice the Golden ratio 
a ~ (75 - l)/2 ^ 0.618 (see Fig. 2). 
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FIG. 2: Power spectral density of the error term in modified Mangoldt function b(n) in comparison 
to the power law 1//^", with the Golden ratio a = {Vb — l)/2. 



The modified Mangoldt function occurs in a natural way from the logarithmic derivative 
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of the quotient 

' n>l 



since —^Jrr = -^i This replaces the relation from the Riemann zeta function where 



Z(s) 

CM = V ^ 



3. THE HYPERBOLIC GEOMETRY OF PHASE NOISE AND 1// FREQUENCY 
NOISE 

The whole theory can be justified by studying the noise in the half plane H = {z = 
u + ly, = —1, y > 0} of coordinates ^ = ^ and y = > ^ and by introducing the 
modular transformations 

The set of images of the filtering line z = u + I under all modular transformations can be 
written as 



k-(- + 7^)l = i- (9) 



q^ 2q!^^ 2qt 

Equation © defines Ford circles jsf [ld| (see Fig. 3) centered at points z = — + ^ with radius 



2^:^ To each 1^ belongs a Ford circle in the upper half plane, which is tangent to the real 
axis at z/ = ^. Ford circles never intersect: they are tangent to each other if and only if they 
belong to fractions which are adiacent in the Farey sequence t < ■ ■ ■ — < Si±pi < 22 ... ^ !_ 

° •J ^ 1 gi 131+92 92 1 

For modular transformations one easily calculates '^'y(z) = — r— nj, ^^^r^ = 7 — r-rw- 

where the symbol means the imaginary part. This implies \^^-^^\{'<s^{z))~'^ = {^z)~^ and 

the invariance of the non- Euclidean metric — = ^'^^ . 

y y 

A basic fact about the modular transformations (jH)) is that they form a discontinuous 
group r SL{2, Z) / {±i }, which is called the modular group. The action of F on the half- 
plane H looks like the one generated by two independent linear translations on the Euclidean 
plane, which is equivalent to tesselate the complex plane C with congruent parallelograms. 
Here one introduces the fundamental domain of F (or modular surface) F = {z & H : \z\ > 

W\ ^ f}; and the family of domains {7(F), 7 G F} induces a tesselation of if jl^. 

Corresponding to the Riemannian metric there is the non-Euclidean Laplacian A = 
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FIG. 3: Ford circles: the mapping of the filtering line under modular transformations 0- The 
arrows indicates that Ford circles were used as an integration path by Rademacher to compute the 
partition function p{n). 
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FIG. 4: The phase angle 6{k) for the scattering of noise waves on the modular surface. 



y'^{-^ + and one can write a non- Euclidean wave equation [l 



ih 



'dt 



— ^(A + -)^. 



Stationary solutions \& = (f>{i^,y) exp(— /cut) satisfy 



A + - + = 0, with P = — — , E = nw. 
4 



(10) 



(11) 



They are special //-independent solutions 4> = y'^ with s = | + Ik, and any other wave can 
be transformed by ((HI) to form the new solution 

y^ 1 



I , T^' s = - + Ik. 



(12) 
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In Gutzwiller's paper j3] there is a detailed discussion on the geometry of solutions 
(fT^ . The i/-independent solution represents a wave propagating in the y-direction = 
y^^"^ exp {Ik In y — lujt). The factor y^^"^ comes from restricting the total flux in a vertical 
strip of constant width such as < < 1, so that with the hyperbolic metric one has 
J^du/y = 1/y. The phase factor is explained by looking at wavefronts y = C^^ and the 
distance proportional to them is dy/y = Iny. 

For equation (fT^ the wavefronts are obtained by assigning to the quantity , — ^-712 some 
constant value that one can choose to be the unity. These wavefronts are circles of equation 
(z/ + + {y — = They are tangent to the real axis at z/q = — ^ and of radius 
Up to a constant shift along the real axis z/, they are the same as the Ford circles shown on 
Fig. 3. 

Thus an outgoing plane wave can be said to start at y = 1, and move up into H with 
increasing y, its wavefronts parallel to the real axis. All other plane waves start at some circle 
which touches the real axis at some arbitrary point z/q of radius Then they contract 
by shrinking their radius while maintaining their point of contact z/q. These wave fronts 
are horocycles in the hyperbolic geometry of the half-plane. They are perpendicular to the 
geodesies, which are Euclidean half-circles, hitting the real axis at uq, at a right angle. 

In addition to (fT^ they are general solutions 13| of the form 

Uz)= ^=l + Ik, (13) 

where Too is the stabilizer in the modular group F of the cusp at infinity, i.e. the subgroup 
of integer translations z ^ z + n. The solution corresponds to waves scattered again and 
again from the modular surface F. For arbitrary s the series fll3|) can be rewritten as 



z 



l\2s 



^y^ ^(Pi qi)=i ip z+q'l^" with the summation performed over all coprime numbers 



iPiili) = 1- The convergence is ensured for 3?(s) > 1. 
The series (tT3|) satisfies a functional equation 



e(2s)0,(z)=e(2-2s)0i_,(2;), (14) 
where ^(s) = 7r"'*/^r(s/2)C(s) is the completed Riemann zeta function. It follows that (fT^ 
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is expanded as 

with S{s) = A{s)Z{s), 
Zis) = '-^, A(.) = MgzlM, (15) 

and the remaining term vanishes exponentially when y — > oo. For s = | + Ik, that is s on 
the critical line, the remainder is Ts{y) = 0. In such a case the scattering coefficient equals 

that is the flux of the reflected wave is equal and opposite sign the incoming flux. The phase 
angle represented on Fig. 4 is a very complicated function of k and is even considered as a 
prototype for quantum chaos jl^ . 

It is observed that the term + Ik) in (fT3j) increases while \Z{^ + Ik)\ decreases 
monotonously to preserve the modulus 1 of S{-^ + Ik). All the stochastic dynamics is 
encoded in the function + Ik). Fig. 4 represents the phase angle 6{k) in the scattering 
coefficient + Ik) = exp{I9{k)) in comparison to the one 6'o(fc) of Z(| + Ik). Apart for 
weak changes the two curves looks similar; this is confirmed from the power spectral density 
which is about the same for the two cases. 

Looking at the logarithmic derivative of the scattering coefficient S{s) = exp{i9{s)), one 
can get the counting function 

dlnS(s) 1 ^, 

e'{k) = — at s = - + Ik. (17) 

which is related to the stochastic factor Z{s) as 

_|M.fM.,r,-.-.B(,),, (18) 

with B{t) = Xl^i^l^)' ^{^) — ^{n)(t){n) /n. By inverting the Mellin transform in (fT^ one 
gets 

Bit) = / -^ -ds. (19) 

^It^ Ju(s)~ioo Z{s) s 

This extends the calculation performed in ^ for getting the error term in the summatory 
Mangoldt function tp{t). There the Riemann zeta function ({s) replaces the quotient Z{s) = 
({s) I C(s + 1) given in ((7j). As reminded in Fig. 2, the error term in the summatory modified 
Mangoldt function is very close to a 1//^" noise, with a the Golden ratio. 
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4. CONCLUDING REMARKS 



From its definition 1// noise is attached to the use of the fast Fourier transform (FFT). 
But the FFT refers to the fast calculation of the discrete Fourier transform (DFT) with 
a finite period g = 2', / a positive integer. In the DFT one starts with all q^^ roots of 
the unity exp{2i7cp/q), p = 1 . . .q and the signal analysis of the arithmetical sequence x{n) 
is performed by projecting onto the n^^ powers (or characters of Z/qZ) with well known 
formulas. 

The signal analysis based on the DFT is not well suited to aperiodic sequences with many 
resonances (by nature a resonance is a primitive root of the unity: {p, q) = 1), and the FFT 
may fail to discover the underlying structure in the spectrum. We recently introduced a new 
method based on Ramanujan sums 3, 3]. 

Cq{n) = exp(2i7r--n,) with {p,q) = 1, (20) 
p=i ^ 

which are n^^ powers of the q^^ primitive roots of the unity. The sums are evaluated from the 
use of Mobius transforms as Cq{n) = ji [-r^) ./'^f \ ■ Here (g, n) is the greatest common 
divisor (g, n) of q and n and fi{n) is the Mobius function. It is if contains a square, 1 if 
n = 1 and (—1)'^ if n is the product of k distinct primes. 



The sums are quasiperiodic versus the time n (ci = 1; C2 = —1, 1; C3 = —1, —1, 2, where 
the bar indicates the period; for example 03(4) = —1) and aperiodic versus the order q of 
the resonance. In particular Cq{n) = fi{q) whenever {n,q) = 1. 

Mobius function can be considered as a coding sequence for prime numbers, as it is the 
case of Mangoldt function. Mangoldt function is related to Mobius function thanks to the 
Ramanujan sums expansion found by Hardy 

tin) = MA(n) = ±tM.i„). (21) 

n ^ 0(g) 

We call such a type of Fourier expansion a Ramanujan- Fourier transform (RFT). General 
formulas are given in our recent publication 3 and in the paper by Gadiyar p] . This author 
also reports on a stimulating conjecture relating the autocorrelation function of b{n) and 
the problem of prime pairs. In the special case (I2H), it is clear that fJ'{q)/4>{q) is the RFT of 
the modified Mangoldt sequence b{n). 
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FIG. 5: Ramanuj an- Fourier transform (RFT) of the error term (upper curve) of modified Mangoldt 
function b{n) in comparison to the function fi{q) / (p{q) (lower curve). 



Using Ramanuj an- Fourier analysis the 1//^" power spectrum gets replaced by a new 
signature shown on Fig. 5, very close to fJ'{q)/(p{q) (up to a scaling factor). There is thus a 
deep relationship between 1// noise in phase locking, the Golden ratio, the Mobius function, 
the modified Mangoldt function, the frequency of windings around the modular surface, and 
the Riemann hypothesis. All these ingredients arise from the hyperbolic geometry of the 
half-plane. 
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